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! Abstract 
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■ 

We show the three-loop integrabihty of large plane-wave matrix theory in a subsector 
Qh! of states comprised of two complex light scalar fields. This is done by diagonalizing the 
^ ■ theory's Hamiltonian in perturbation theory and taking the large limit. At one- loop level 
the result is known to be equal to the Heisenberg spin-1/2 chain, which is a well-known 
^ I integrable system. Here, integrability implies the existence of hidden conserved charges and 
5^ ! results in a degeneracy of parity pairs in the spectrum. In order to confirm integrability 
at higher loops, we show that this degeneracy is not lifted and that (corrected) conserved 
charges exist. Plane- wave matrix theory is intricately connected to A/" = 4 Super Yang- 
Mills, as it arises as a consistent reduction of the gauge theory on a three-sphere. We find 
that after appropriately renormalizing the mass parameter of the plane-wave matrix theory 
the effective Hamiltonian is identical to the dilatation operator of A/" = 4 Super Yang-Mills 
theory in the considered subsector. Our results therefore represent a strong support for the 
conjectured three-loop integrability of planar A/" = 4 SYM and are in disagreement with a 
recent dual string theory finding. Finally, we study the stability of the large N integrability 
against nonsupersymmetric deformations of the model. 



1 Introduction and Conclusions 



Recently a number of fascinating developments have pointed to the existence of integrable 
structures in the large limit of four- dimensional A/" = 4 supersymmetric SU(A^) Yang- 
Mills theory |ll|2llSlll] . These developments have been paralleled by investigations on the 
integrability of the dual AdS^ x string [S]. In this paper we shall turn to another setting 
of string/M-theory in which integrable structures emerge: The large N limit of plane-wave 
matrix theory [Hj. This SU(A^) matrix quantum mechanics represents a mass deformed 
version of the "BFSS" matrix theory 7J, thereby preserving all of its supersymmetries. In its 
large limit it is believed to describe the microscopic degrees of freedom of M-theory in the 
maximally supersymmetric plane- wave background jHj. Opposed to its fiat space cousin the 
model enjoys a number of pleasant features: It possesses a discrete spectrum, the model may 
be treated perturbatively in the large mass limit P , and its spectrum contains an infinite set 
of exactly protected states due to its underlying SU(2|4) symmetry structure [1211111112] • 
The present worlc, however, is mainly motivated by the direct connection of plane-wave 
matrix theory to A/" = 4 Super Yang-Mills on R x 5^. The matrix quantum mechanics 
arises in a consistent reduction of the four-dimensional superconformal gauge theory on 
the • As the scaling dimensions of operators in a conformal theory on are related 

to the energy spectrum on R x S"^ in radial quantization, the dilatation operator of the 
four-dimensional Super Yang-Mills theory should bear a relation to the Hamiltonian of 
the plane-wave matrix model. And indeed in ^Hl it was shown that at the one-loop level 
the spectrum of anomalous dimensions of gauge theory operators agrees with the leading 
perturbative corrections to the spectrum of associated states in the plane-wave matrix 
quantum mechanics. An immediate consequence of this observation is, that the one-loop 
integrability of the Super Yang-Mills model, to be discussed below, carries over to plane- 
wave matrix theory. 

Integrability in A/" = 4 Super Yang-Mills shows up in the perturbative study of the di- 
latation operator, whose eigenvalues yield the spectrum of scaling dimensions of composite 
operators. In principle the scaling dimensions are extracted from diagonal two-point func- 
tions, however, in jElHj a new efficient technique was introduced to directly extract the 
dilatation operator order by order in a loop expansion from the gauge theory. This led to 
the remarkable discovery of Minahan and Zarembo Hf] of an exact equivalence between the 
planar, one-loop piece of the dilatation operator (in the pure scalar excitation sector) with 
the Hamiltonian of an integrable 5*0(6) spin chain. In this picture the spins on the chain 
are given by the scalar fields (pi (z = 1, . . . , 6) of A/" = 4 Super Yang-Mills thread together 
in a single trace structure, tr(0j^ . . . 0j„). This finding was subsequently generalized to the 
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full Super Yang-Mills field content |3j establishing an equivalence of the planar one-loop 
dilatation operator with the Hamiltonian of an integrable SU(2,2|4) spin chain. 

Obviously one would like to know whether this integrability is a mere artifact of the 
one-loop approximation, or whether it represents a genuine hidden symmetry of planar Su- 
per Yang-Mills. To answer this question higher loop considerations suggest themselves. For 
this it is advantageous to restrict one's attention to a minimal subsector of field excitations, 
which does not mix in perturbation theory with other excitations. Such a sector is given by 
the SU(2) excitations of two complex scalars, which we denote Z and 0. The correspond- 
ing one-loop planar dilatation operator is then identical to the famous XXX1/2 Heisenberg 
quantum spin chain model, H = ^^=1 Sj ■ Sj^i, where Sj is the spin 1/2 operator acting 
on the jth site of the chain, being in an up (Z) or down (0) state. This model of a one di- 
mensional magnet is known to be integrable, i.e. one may construct a set of L — 1 mutually 
commuting charges which contain H ^Hj- In [2] a two-loop computation in this SU(2) sec- 
tor of A/" = 4 Super Yang-Mills was performed showing that integrability is preserved. The 
higher loop contributions to the planar dilatation operator result in deformations of the 
spin chain Hamiltonian incorporating spin interactions beyond nearest neighbours, i.e. at / 
loop level one is dealing with the interactions of / + 1 spins. The integrability of the higher 
loop deformation manifests itself in the preservation of degeneracies in the spectrum, which 
may be attributed at one-loop level to the hidden symmetry. However, a generalization of 
the powerful methods known from integrable spin chains for the construction of commuting 
charges at higher loop orders remains unknown, nevertheless the first few charges beyond 
the deformed Hamiltonian were constructed explicitly. It is important to note that while the 
one-loop dilatation operator is exactly integrable, the combined one- and two-loop dilata- 
tion operator is only integrable up to terms of three-loop order. Hence only the full, planar, 
all-loop dilatation operator can display exact integrability - a fascinating perspective on 
the large N spectrum of A/" = 4 Super Yang- Mills. 

Based on the assumption of three-loop integrability the authors of [2] where actually 
able to determine the form of the planar three-loop dilatation operator. However, this con- 
jecture was recently challenged by the work of Callan, Lee, McLoughlin, Schwarz, Swanson 
and Wu JH] from the dual string theory side. These authors computed the corrections to 
the spectrum of IIB plane-wave superstring theory due to the presence of first curvature 
corrections that arise in the expansion of AdS^ x about the plane- wave limit. These cor- 
rections translate into certain 1/J corrections of the scaling dimensions of BMN operators 
in the BMN limit jH] and can be compared to the two-loop result and three-loop conjecture 
of Beisert, Kristjansen and Staudacher Indeed the two-loop corrections of the gauge 
theory match with the string result, however, at three-loop order the conjectured gauge 
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theory result disagrees with the string based findings. Hence, if the involved analysis of ^H] 
is indeed correct, then string theory predicts the breakdown of integrability in A/" = 4 Super 
Yang-Mills at the three-loop level! To clarify this issue from a gauge theory perspective a 
direct three-loop analysis would be very desirable. Fortunately the aforementioned relation 
of A/" = 4 Super Yang-Mills to plane-wave matrix theory enables one to perform such a 
three-loop computation in a simpler quantum mechanical setting. This is the goal of our 
paper. 

Our analysis consists of an essentially straightforward application of degenerate quan- 
tum mechanical perturbation theory to the third non-vanishing order, which is equivalent 
to a field theoretic three-loop computation. Needless to say that this has been only pos- 
sible by making massive use of symbolic computer algebra systems. The outcome of our 
analysis is positive for integrability: Plane-wave matrix theory remains integrable up to 
the three-loop level in the minimal SU(2) sector mentioned above. Technically this means 
that the degeneracy of "planar parity pairs" noted in extends to the three-loop order. 
Also the first charge beyond the Hamiltonian can be explicitly constructed. Clearly this 
does not prove the integrability of the full four- dimensional theory, as we are dealing with 
a reduced model here, but it does make it very likely. This is furthermore supported by the 
following remarkable observation. As found in [131 the agreement of the Super Yang-Mills 
dilatation operator with the perturbative effective plane-wave matrix theory Hamiltonian 
ceases to exist at two-loop order. However, if one allows for a renormalization of the rela- 
tion between the (dimensionless) mass parameter M of matrix theory and the Yang-Mills 
coupling constant gyu of the form^ 

M3 327r2 V 32 7r2^^^^ 252 tt^ ^™ ^ ^^^^ V ' 

then this agreement extends to three-loop order! That is our three-loop matrix theory com- 
putation reproduces the conjectured three-loop dilatation operator of Beisert, Kristjansen 
and Staudacher [2j with the above choice of renormalized mass parameter. Indeed such a 
renormalization is to be expected if one were to truly path integrate out the higher Kaluza- 
Klein modes of the gauge theory on to obtain an effective one dimensional theory of the 
zero-modes. The resulting effective theory would then be given by the plane-wave matrix 
model plus higher loop corrections in A = g^-^ N and 1/A^. These quantum loop corrections 
should manifest themselves in a renormalization of the mass parameter M of the plane- 
wave matrix model and - but our findings indicate that this is not the case - the addition 
-^This observation grew out of discussions with N. Beisert and M. Staudacher. 
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of higher order interaction terms. It would be very interesting to study this in more detail, 
in particular the role of (maximal) supersymmetry in this non-emergence of higher order 
interaction terms. 

In addition to that, we have explored how stable the integrability is to deformations of 
the theory. It turns out that supersymmetry is absolutely unessential for integrability in 
the matrix model up to three-loops. Switching off all fermions in the plane-wave matrix 
theory does not destroy integrability. Furthermore deformations of the bosonic field content 
are also possible without loosing integrability, as long as the structure of Yang-Mills type 
interactions is maintained. Our findings suggest that the hidden integrability, if it exists 
to all loops, is a genuine property of large N gauge quantum mechanics. It would be very 
nice to prove integrability in these models. To this end a formalism enabling one to directly 
work in the planar limit of the theory is needed. Maybe the structures introduced in ^Jj 
could be of use here. 

Naturally a further important question is to understand how these (potential) integrable 
structures are to be interpreted in M-theory. 

Finally, let us mention that an alternative method to address the question of three-loop 
integrability of the full four-dimensional TV = 4 gauge theory lies in the study of how 
strongly the underlying SU(2, 2|4) superalgebra - or a suitable subgroup thereof - restricts 
the three-loop structure of the dilatation operator. This will be addressed in jTHj. 



2 Plane-wave Matrix Theory 

The Hamiltonian of plane-wave matrix theory can be written as if = Hq + V1 + V2 with 



Hn = tl 



\PiPi + \ (f )' X,X, + \M^X,X, - ^z^7i23^ , (2a) 



Vi = -Mieej'k' tr Xi,Xj,Xk> - tr 07/[X/, 9] , (2b) 
V2 = -^tT[Xj,Xj][Xj,Xj] . (2c) 

The degrees of freedom are the light scalars Xj {i = 1, ... ,6) of mass M/2, the heavy 
scalars Xj/ {i' = 1, . . . , 3) of mass M and a SO (9) Majorana spinor 6a {en = 1, ... , 16) of 
mass 3M/4 ^. All fields are N x N traceless Hermitian matrices by virtue of the gauge 
group SU(X). The index / = embraces all scalars. In the fermionic sector we work 

^The mass parameter here is a dimensionless quantity given by M = in terms of the M-theoretic 
quantities /i (mass parameter of the plane-wave background, see eq. (1) of |13p. Ip (eleven dimensional 
Planck length) and R (radius of compact eleventh dimension). 
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in a representation with charge conjugation matrix equal to unity and with symmetric 
euchdean Dirac matrices (7/)a/3. Note that it is the bosonic Xj sector which is intimately 
related to the scalar field sector of A/" = 4 supersymmetric Yang- Mills theory. 
We introduce the following creation and annihilation operators 



f iai + a] 



Xji 



/2M 



h, - h], 



with 6^ := , 



where 11^ := |(1 ± i'jus)- They satisfy the canonical ( ant i-) commutation relations 



(3) 



(4) 



and the free Hamiltonian now reads 



Ho = tT 



—ajai + Mb],bi' + —OtO- 



(5) 



Physical states are constraint to be gauge invariant and are given by traces over words in 
the creation operators. 



3 Perturbation Theory 

We are faced with an eigenvalue problem of the form 

{Ho + kVi + «:V2)|0(/€)> = E{k)\<P{k)), (6) 

where the spectrum of the free (k = 0) Hamiltonian Hq is known. Let us sketch how (jH)) 
is solved perturbatively in k and how we deal with potential degeneracies. Say we are 
interested in the perturbative shift of a specific L-fold degenerate free energy value Eq. 
Then we choose an arbitrary basis i = 1, . . . , L} for all states with this free energy 

i/o|0.> = ^o|0.> (z = l,...,L). (7) 

For later convenience we do not assume this basis to be orthonormal. All other states are 
labeled by In) and their free energies by En. These energies may, of course, be degenerate 
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as well. Here, however, we define the states in such a way that 

{n\m) = 6nm , {n\(pi) = , (8) 

which is always possible. 

When the interaction is turned on (k ^ 0), the states get corrected to 
Clearly, the corrected states [^^(k)) are generically not eigenstates of the full Hamiltonian 
H = Hq + kVi + K^V2 since there may occur mixing among the states However, we do 
not want to deal with mixing problems and the possible lift of the original degeneracy at 
this stage, and therefore introduce an energy matrix Eij{K,): 

{Ho + kVi + /€V2)|0,(/€)> = Ei,iK)\(l),{K)) . (9) 

Later, this matrix can be diagonalized and its eigenvalues are the final corrected energies. 
In this way we have decoupled the mixing of states in the degenerate subspace from the 
admixture of states from outside this subspace, which turns out to be extremely convenient 
for practical computations. 

Now we wish to determine an energy operator T(k) whose matrix elements give exactly 
the energy matrix: 

Eijin) = {4>,\TiK)\(Pi) . (10) 
In |TU|) we have introduced the dual basis i = 1, . . . , L} defined by 

= '^u- • (11) 

This is necessary because we did not demand the orthonormality of On the other 

hand it is also very convenient for the following reason. We can start from any arbitrary 
basis and do not need to worry about orthogonalizing it but may immediately apply 

T(k) with a result of the form 

T(/s:)|0,) = + J2^^nin)\n) . (12) 

j n 

From the expansion coefficients we then simply read off the energy matrix as 

E,,{K) = U,{n) . (13) 
In order to find T{k) we follow a standard procedure. States and energy matrix are 
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expanded in k: 



|0.(«:)> = 10.) + K (4^V.-> + h>) + (4'V.-> + h>) + • • • (14) 
E,,{^) = Elf + KE\f + K'Ef + ... with Ei;) = (15) 

Plugging these expansions into Q yields one equation for each power of k. When the 
resulting equations are projected onto and in\ one obtains expressions for E\'^'^ and 
if)^^ . The coefficients -^^ remain undetermined and may be set to zero. Iteratively we then 
find the energy operator 

oo 

T{n) = Y,^''Tk (16) 



k=Q 



up to third loop order 



To = Ho (17a) 

T2 = Vi^Vi + V2 (17b) 

T4 = ViA^iAViAri + ^lAViAVa + l/iA^aAVi + l/sAV^iAVi + l/aAKa (17c) 
- ViAViPTs 



Te = V^iAViA\/iA\/iAViAl^i 

+ 14A\/iA\/iAViAy2 + l^iAVlA\/iAV2AVi + l/iA^iAyaAFiAK 

+ I/1AV2AV1AV1AV1 + l/aA^iAViAriAri 

+ ■l/iAyiA\/2AV2 + 14Ay2AViA\/2 + VlAy2A\/2A\/i 

+ 1/2AViAVlA1^2 + V'2AViAV2AVi + ■l/2AV^2AViAVl 

+ V2^V2^V2 

- (yiAViA\/iAFi + FiA\/iAViAFi + yiA^iAViAVi 
+ I/1AV1AV2 + ViAV2Ayi + I^AViA^i 

+ l^iAyiAV2 + ViAV2^^Vi + F2AVIAV1 

+ V2AV2)PT2 

- ViA^ViPTi + ViA^ViPT2PT2 (17d) 
giving rise to the expansion of the energy matrix Moreover vanishes for odd /c. 
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Here we have introduced the "propagator" and the projector 



A = T^ TT and P = Y}<^>^){^. (18) 

respectively, where P projects onto the subspace of states with free energy Eq. 



4 Energy Operator in SU(2) subsector 

In this section we compute the energy operator T up to third loop order for the SU(2) 
subsector of the matrix model which is spanned by gauge invariant states composed of the 
two complex fields 

Z=—={a\ + ia2) and (j) = —={a\ + ia}^ . (19) 

V 2 V 2 

The arguments given in section 3.1. of P] for the group theoretical constraints on mixing of 
the analogue SU(2) sector in A/" = 4 supersymmetric Yang- Mills field theory are inherited 
by the plane-wave matrix theory. This is so as these arguments are merely based on the 
representation content of the excitations with respect to the SO (6) /^-symmetry group 
and the classical scaling dimensions of the involved fields, which directly carry over to the 
matrix model. In terms of the SO (6) Dynkin labels the states we consider transform in the 
representations [p, g, p] and have free energy Eq = ^ (2p + q) . 

We apply a procedure similar to our previous two-loop computation jTHj. We normal 
order the terms given in ()17p and discard all pieces which annihilate multi-trace states of 
Z and (p. Furthermore we can neglect all terms which are not relevant for the expansion 
coefficients tij but only for tin in (|T^ . The relevant terms are exactly those which conserve 
the free energy, i.e. the number of fields of a state. 

In addition to the above fields, the energy operator consists also of the canonically 
conjugated fields 

Z = —^{ai — ia2) and (p = —i={a3 — ia^) , (20) 

V 2 V 2 



which satisfy 



[Zrs, Ztu] — [4'rs,4'tu] — ^st^ru ~ ^^rs^tu ■ (21) 



Since our intention is to compare the energy operator of the matrix theory to the dilatation 
generator of the field theory, we make the following definitions. Perturbation theory is done 
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for large masses in an expansion in As shown in |13j, this mass parameter should be 

related to the SYM coupling constant Qym by 



M3 327r2 V 2 2N 

in order to match all one-loop results in the scalar S0(6) sector of both theories. In ()22|) we 

have also given the relation to two further shorthands G := and A := G'^N. Now we 

An 

define the counterpart of the dilatation generator from the energy operator T by factoring 
out the free energy M/2 of an SO (6) scalar, i.e. the mass of Z and 0: 

M M °° 

T{A)=:-D{A)=:-J2^'D,,, (23) 

k=0 

which explicitly means 

2 A'P 
^o = ^To, D, = —T,, ^4 = ^T,, D, = ^n. (24) 

We will refer to -D(A) as the dilatation operator and to its matrix elements (0j|Z^(A)|(/)j) 
as the dilatation matrix. The eigenvalues of -D(A) will be denoted as A (A). 

Performing the computation using Form ^H] as well as Mathematica, we find the fol- 
lowing result: 

Do = : tT{ZZ + 00) : , (25a) 
D2= - ^:tr[Z,0][Z,0]:, (25b) 

D4=^:tr[Z,0][Z,0]: 

+ J^: (tr[Z,0][Z,[Z,[Z,0]]]+tr[Z,0][0,[0,[Z,0]]]) :, (25c) 
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115 



N 



:tT[Z,4>][Z,^] : 



+ 



61 
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: {trm, Z],T%Z][Z, [T", [0, Z]]] + tr[[[0, Z], T"], 0] [0, [T", [0, Z] 



16iV3L(^^^[^"' 



+ 



80 
31 



195 



16A^3 
111 



: (trZ[T^ 



31 :(tr^[T^ 
+ f : (trZ[^^ 

-j3^M0,^][[ 



+ tr[[[0,Z] 



-tr[[[0,Z] 



[z, [T^ [0, [0, z]]]]]] + tr 0[^^ [0, [0, [^^ [z, [z, 

Z, [0, [T^ [Z, [0, Z]]]]]] + tr 0[T^ [0, [Z, [T^ [0, [Z, 
0, [Z, [T^ [Z, [0, Z]]]]]] + tr 0[T^ [Z, [0, [T^ [0, [Z, 

z, [z, [^^ [0, [0, z]]]]]] + tr 0[^^ [0, [0, [T^ [z, [z, 

Z, [0, [T^ [Z, [0, Z]]]]]] + tr 0[T^ [0, [Z, [T", [0, [Z, 

0, [z, [^^ [z, [0, z]]]]]] + tr 0[T^ [z, [0, [^^ [0, [z, 



Z,0],0],[Z, [Z,0]]]: 



Z], Z] [Z, [Z, [Z, 0]]] + tr[[[0, Z], 0], 0] [0, [0, [Z, 
0], 0][Z, [Z, [Z, 0]]] + tr[[[0, Z], Z], Z][0, [0, [Z, 0]]]) : 
Z], 0] [0, [Z, [Z, 0]]] + tr[[[0, Z], 0], Z] [Z, [0, [Z, 
Z], Z][0, [0, [Z, 0]]] - tr[[[0, Z],Z], 0][Z, [0, [Z, , 



5 Planar limit 



(25d) 



We now take the planar limit of the above SU(2) dilatation operator D, i.e. we derive a 
new operator (again called D), which yields the dilatation matrix in the limit N —* oo with 
A = fixed. Since in the planar limit single and multi-trace states do not mix, we may and 
we do restrict the action of D to the closed subset of single trace states. 

Hence, the dilatation operator in the planar limit may be written as a sum of permu- 
tations of the fields in the trace. We use the notation of Let Pk^M exchange the fields 
at sites ki and k2 in the trace (where the sites are periodically identified) and define 



{ni,n2,...} := ^ Pk+ni,k+ni+lPt 



k+n2,k+n2+l ' ' ' t 



(26) 



k=l 



where L is the number of fields the operator is applied to, i.e. the length of the trace 
respectively chain. The identity operator is denoted by {}. There are many trivial relations 
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such as 



{...,n,n,. ..} = {.. .,...} (27) 
{... ,n,m, ...} = {... ,171,11, .. .} for \n — m\ > 2 (28) 

{111,112, ...} = {ill + m, 112 + m, .. .} (29) 

as well as the identity 

{...,...} + {...,11 ±l,n,...} + {...,11,11 ±1,...} 

-{...,11,...} -{..., n±l,... }-{..., n,n±l,n,...} = (30) 

which only holds in the SU(2) subsector and expresses the fact that two objects cannot be 
placed completely antisymmetric onto three sites. 

In order to find the planar version of the dilatation operator, we have made an ansatz 
for all independent permutations and determined their coefficients such that the N ^ oo 
limit of ()25p is correctly captured. We find the rather compact result: 

^0 = {} , (31a) 

D2 = 2{}-2{0}, (31b) 

D, = -15{} + 19{0} - 2({0, 1} + {1, 0}) , (31c) 
Dq = 187{} - 259{0} + 38({0, 1} + {1, 0}) 

+ 4{0, 2} - 4({0, 1, 2} + {2, 1, 0}) - 2({0, 2, 1} - {1, 0, 2}) . (31d) 

The presence of the antisymmetric term ({0, 2, 1} — {1, 0, 2}) in the three-loop part might 
appear unexpected. However, it is just a consequence of the fact that the energy operator 
()17|) is not Hermitian^. By a change of basis one can remove this term: after the following 
similarity transformation 

D'{A) := e-^^°>D(A)e^^°> (32) 
we obtain D'q = D'q, D'2 = D2, D'^ = D4 and 

D', = 187{} - 259{0} + 38({0, 1} + {1, 0}) + 4{0, 2} - 4({0, 1, 2} + {2, 1, 0}) . (33) 

In [13] we observed that the two-loop dilatation operator ()31c|l failed to agree with the 
two-loop dilatation operator of SYM if one relates M to Qym according to (j221)- However, 



^As a matter of fact T is quasi-Hermitian, i.e. there exists an operator S such that = S S. 
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if one allows a renormalization of this relation, then agreement with the full (planar and 
non-planar) two-loop dilatation operator of SYM can be achieved. Moreover, we discover 
that this is true for the planar three-loop operator as well. Let thus Gr be the renormalized 
coupling related to G by 



G^ = Gl + —Gl - WN'^Gl + 0{Gl) . (34) 
Then the dilatation operator D' written in terms of := Gf.N becomes 

D'r,o = {} > (35a) 

D;,2 = 2{} - 2{0} , (35b) 

^;,4 = -8{} + 12{0} - 2({0, 1} + {1, 0}) , (35c) 



D'rfi = 60{} - 104{0} + 24({0, 1} + {1, 0}) + 4{0, 2} - 4({0, 1, 2} + {2, 1, 0}) , (35d) 

which is exactly the proposed field theory result, cf. equation (F.3) of P). There the planar 
three-loop dilatation operator of SYM was determined using integrability as an input. Here 
we find it without any assumption. 

6 Integrability 

Since the planar dilatation operator (jHT|) of plane-wave matrix theory in the SU(2) subsector 
can be transformed to the corresponding dilatation operator of SYM, we know from the 
investigations of j.2, that it is integrable. I.e. there exists a charge 

oo 

U{A) = J2^'U2k (36) 

fe=i 

which commutes with the dilatation operator. This charge pairs up states of different parity 
and is therefore responsible for their degenerate energies. Here parity refers to the operation 
of inverting the order of fields in a trace, e.g. P ti {Z^(f)^Z(f)) = tT{(j)Z(j)^Z^). The dilatation 
operator commutes with parity, which is the reason for states of opposite parity not mixing. 
A priori there is no reason for a relation between the spectra of positive and negative parity 
states. In the matrix model, however, we discover that there is a maximal degeneracy of 
these spectra which we have explicitly evaluated up to dimension 10: Each representation 
that admits both parities has maximal degeneracy, see figure 1. 

It is instructive to see this explicitly in the two lowest dimensional representations with 
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[5,0,5] 

[4,0,4] [4,1,4] [4,2,4] 

[3,0,3] [3,1,3] [3,2,3] [3,3,3] [3,4,3] 

[2,0,2] [2,1,2] [2,2,2] [2,3,2] [2,4,2] [2,5,2] [2,6,2] 



4 5 6 7 8 9 10 

Dimension A(A = 0.125) 

Figure 1: The complete spectrum of unprotected states in the SU(2) subsector up to bare dimen- 
sion 10 at three-loop order for A = 0.125. States of positive (negative) parity are depicted as dark 
(light) lines. The states are grouped according to the Dynkin labels of their representation. Note 
that whenever states of positive and negative parity exist within a representation, then they are 
degenerate. The explicit results for the energy shifts are spelled out in the apppendix. 

S0(6) Dynkin labels [3, 1,3] and [3,2,3]. In the [3, 1,3] representation of Aq = 7 one has 
two single-trace operators of opposite parity 

^[3,1,3] ^ 2tr[000ZZZZ] - 3 tr[00Z(^ZZZ] + 2 tr[00ZZ0ZZ] 

- 3tr[00ZZZ0Z] + 2tr[0Z0Z0ZZ] , 
^[3,1,3] ^ _ tj.[00Z0ZZZ] + ix[(t)(t)ZZZ(t)Z] . 

The associated energy shifts up to order read 

5A+ = 5A_ = lOA^ - 65A^ + 765A^ (37) 

and are degenerate. 

Similarly for the representation [3, 2, 3] of Aq = 8 one finds one single-trace operator of 
positive parity 

^[3,2,3] ^ i^[(i)(f)(i)ZZZZZ] - ix[(t)(t)Z(t)ZZZZ] - tr[(l)(f)ZZZZ(l)Z] + tr[(l)Z(l)ZZ(f)ZZ] 
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and two operators of negative parity 



^[3,2,3] ^ -tr[(f)(t)Z(t)ZZZZ] + tr[(t)(t)ZZ(t)ZZZ] - tr[00ZZZ0ZZ] + tr[00ZZZZ0Z] , 
^[3,2,3] ^ _ tr[(/,0Z0ZZZZ] - ti[(l)(l)ZZ(l)ZZZ] + tr[</)0ZZZ0ZZ] + tr[0</)ZZZZ0Z] . 

In order to compute the energy shifts in the negative parity sector one has to diagonahze 
a 2 X 2 matrix, C_ := {0^^f^\0^^f^^y 

D,0^ = °] 0_ B,0_ = ] 0_ = h " V- ■ (38) 

The outcome of this straightforward exercise are the energy shifts 

5A+ = SA^ - 48A^ + 536A*^ (39) 

f SA^ - 48A^ + 536A^ , , 

[12A2-78A4 + 942A6 ^ 

which again display degeneracy of opposite parity states. 

Responsible for this degeneracy is the charge U of which commutes with H. Up 
to three-loop orders it is given by 

t/2 = 4({l,0}-{0,l}), (41a) 
f/4 = 8({2,l,0}-{0,l,2}), (41b) 
f/e = 8({1, 0, 3} - {0, 1, 3} + {0, 3, 2} - {0, 2, 3}) + 18({3, 2, 1, 0} - {0, 1, 2, 3}) 
+ 38{0, 1, 3, 2} - 30{2, 1, 0, 3} - 38{0, 3, 2, 1} + 30{1, 0, 2, 3} 
-34({1,0,3,2}-{0,2,1,3}). (41c) 



7 Nonsupersymmetric Deformations 

In this final section we turn to the question of how strongly the observed three-loop integra- 
bility depends on the details of the model. A number of deformations are easily implemented 
by switching of fields in our computer programs. 

To begin with let us consider the purely bosonic SO (3) x SO (6) model by setting the 
fermions to zero in We calculate the effective vertex for this bosonic model similar to 
the above computation but subtract off the vacuum energy shift - which due to lack of 
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supersymmetry does not vanish - by hand^. Then the planar dilatation generator reads: 



Do^{}, (42a) 

D2 = 14{} - 2{0} , (42b) 

= -558{} + 122{0} - 4({0, 1} + {1, 0}) , (42c) 

D, = 24|29{} - 63|83|0} + 584({0, 1} + {1, 0}) 

+ 16{0, 2} - 16({0, 1, 2} + {2, 1, 0}) - 8({0, 2, 1} - {1, 0, 2}) . (42d) 



It is important to note that the planar parity pairs are degenerate under the individual 
action of the operators 

{}, {0}, {0,1} + {1,0} and 

2(0, 2} - 2({0, 1, 2} + (2, 1, 0}) - ({0, 2, 1} - {1, 0, 2}) . 

Hence the purely bosonic SO (3) x SO (6) model displays integrability up to three-loops as 
well! 

It turns out to be rather hard to break the hidden integrability structure of a large 
matrix quantum mechanical model. For this we have studied a "brutal" deformation and 
considered the bosonic SO{d) model {i = 1, . . . ,d) 

H = —ii a\a, - - tr X^X^X^X^ + ^ tr X,X,X^X^ (43) 

which for parameters a ^ 13 has lost all its traces of a higher dimensional gauge field theory 
origin. Then one finds (again dropping the vacuum energy shift) 

L'o = {} , (44a) 

L>2 = (-4a + (2d + 3)/3){}-2a{0}, (44b) 
Di = \ ((-100 - I2d)a^ + (168 + 96d)a/? + (-59 - 7Ad - 2Ad^)(3^) {} 

- 2 (I6a2 _ (^17 + M)a(5) {0} - Aa^{{Q, 1} + {1, 0}) , (44c) 

*The vacuum energy shift is actually dominant for large A''. 
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De = i ((-560c/ - 2440)a^ + {280(f + 3668rf + 60AA)a^P 

+ (-1536^2 - 5274d - 4592)a/3^ + (256d^ + 1199rf^ + IQlOd + 1029)/5^) {} 
+ ((-36c/ - 832)6^^ + (768c/ + 1700)6^^/5 + {-192d^ - 797c/ - 791)ci/?2) {0} 
+ (-1646^3 + 4(20c/ + 37)a'^P) ({0, 1} + {1, 0}) 

+ 166^2/5(0, 2} - 16a3({0, 1, 2} + {2, 1, 0}) - 8a%{0, 2, 1} - {1, 0, 2}) . (44d) 

We see that the integrabihty of the model ()43|1 only breaks down at the three-loop order, 
as the last line of ()44d|) destroys integrabihty for a j3. If one chooses a = (3, however, 
integrabihty is stable up to three-loops. 

We interpret the findings of this section as an indication that integrabihty might be a 
generic effect in large N gauge quantum mechanical models and independent of supersym- 
metry. In the same instance we learn from the example ()43p . that if integrabihty is to break 
down in a concrete model, it does so at rather high loop-orders (here three). This should 
caution us in our expectancy of an exact integrabihty of plane- wave matrix theory or A/" = 4 
Super Yang-Mills based on three-loop results, even if the idea seems too appealing to be 
false. 
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Appendix 



An 






^A 


4 


[2,0,2] 


izA — yu7v + iuyo/v 




r 




To 1 ol 

[2,1,2] 


oA — ozA + oooA 




u 


\0 9 91 

[3,0,3] 


14.472A2 - 107.01A^ + 1327.4A6 


12A2 - 78A4 + 930A6 


7 


[2, 3, 2] 


4A — 20A + 202. 5A° 






12 A — 84A + 997. 5 A° 






To 1 ol 

[3, 1,3] 


iUA — ooA + 7o5A 


iUA — DoA + 7doA 


8 


[2,4,2] 


3.0121A2 - 13.863A^ + 137.83A6 
15.208A2 - 112.68A^ + 1412.8A6 

Q7Sn9A2 fi'i /I A 4 _i_ 71 Q /1 1 A 6 

y. ( ouziv — oo.4ooiv -|- ( iy.4ii\. 






[3,2,3] 


8A — 48A + 53dA 


8A — 48A + o3dA 






12A2 - 78A4 + 942A6 




[4, 0, 4] 


0.491/A — 3U.275A + 293. 55A 
22.604A2 - 167.80A4 + 2059.3A6 
10.904A2 - 69.927A4 + 859.19A6 
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[2,5,2] 


2.3432A — 10.159A + 100.35A° 






oA2 a o \ 4- 1 rroAfS 

8A — 48A + 552A° 

13.657A2 - 97.841A^ + 1189. 7A6 






[6, 6, 6\ 


fi 4p;'?9A2 Qp; 770 a 1 -4- '^89 1 1 A^ 

11.040A2 - 72.143A^ + 858.37A6 
16.507A2 - 113.09A^ + 1356.5A6 


fi '^^ 77nA'l -1- '^89 1 1 A^ 

11.040A2 - 72.143A'^ + 858.37A6 
16.507A2 - 113.09A^ + 1356.5A6 




[4,1,4] 


5.072A2 - 22.574A^ + 220.97A6 






lOA^ - 65A^ + 785A6 

18.928A2 - 133.43A^ + 1591.0A6 


10A2 - 65A4 + 785A6 



Table 1: Energy shifts for all unprotected states in SU(2) subsector up to dimension Aq = 9. 
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